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ABSTRACT 

Resent observations of a number of galaxy clusters using the Sunyaev-Zel'dovich effect indicate that 
about 1/3 of baryonic mass is missing from the hot intracluster medium (ICM), which is significantly 
larger than the fraction of stars and cool gas, which account for only about 10%. Here we address the 
question whether the remaining 22 ± 10% can be accounted for by thermal evaporation of gas from 
clusters. We have found that evaporation can occur only from the cluster "surface", r ~ rvir, and not 
from it's interior. We evaluated particle diffusion through the magnetized ICM for several scenarios of 
ISM turbulence and found that diffusivity is suppressed by at least a factor of 100 or more, compared 
to the Spitzer value. Thus, only particles from radii r ^ 0.9rvir can evaporate. Diffusion of particles 
from inside the cluster, r <^ O.Qrvir, takes longer than the Hubble time. This lowers the cluster-averaged 
fraction of the evaporated hot gas to few percent or less. However, if the missing hot component is 
indeed due to evaporation, this strongly constrains the magnetic field structure in the cluster envelope, 
namely either (i) the gas is completely unmagnetized {B < 10~^^ gauss) in the cluster halo or (ii) the 
magnetic fields in the ICM are rather homogeneous and non-turbulent. 

Subject headings: galaxies: clusters: general — diffusion — magnetic fields 



L INTRODUCTION 

Observations of the Sunyaev-Zel'dovich effect in galaxy 
clusters directly probe the energy content of the intraclus- 
ter medium (ICM) in the model-independe nt way, and it is 
complementary to X-ray observations (see IVoit |2005| and 
references therein). The study of 193 clusters with tem- 
peratures above 3 keV using the cosmic microvawe back- 
ground dat a from WMAP y ielded a remarkable and puz- 
zling result (jAfshordi. et al.| [20G6). It has been shown that 
SZ effect would account for a half of the thermal energy 
of the ICM provided all the baryons are in the hot phase. 
This result implies that a significant fraction of baryons, 
about 32 ± 10% by mass, is rn issing from the ICM phase. 

It has been proposed fcoebi 120071 ) that gas evaporation 
can reduce internal energy of gas in a cluster. Evaporation 
involves particles from the thermal tail of the Maxwellian 
distribution because their velocities exceed the escape ve- 
locity from a cluster, hence they are gravitationally un- 
bound and can leave the cluster. Since these particles 
are suprathermal, i.e., their speed exceeds the thermal 
speed, the overall energy content of the cluster will be 
decreasing with time. Because the ICM gas is mildly 
collisional, the transport of these particles will be dif- 
fusive. However, in the presence of magnetic fields in 
the ICM, particle diffusion is often substantially reduced 



(iRechester fc RosenbluthlllQTS ^'Chandran fc CowleyllQQS : 



Naravan fc MedvedevI l2001l: Malvshkin fc Kulsrudil2001 



Malvshkinll200ll : Lazarian 2006ll . The reduction coefficient 
fB<^ has been introduced; its value has not been accu- 
rately evaluated, however. 

Magnetic fields in galaxy clusters are highly turbulent 
with the fluctuating component being of order the mean 
field, SB ^ B, and is tangled on small scales. The 
magnetic field spectrum deduced from the rotation mea- 
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sure observations is peaked at about / b ^ one — few kpc, 
which sets the characteristic correlation length of the field 
(|Vogt fc Enfih n 2003, 2005: lEnfilin. et a]||2005[ ). We eval- 
uate the diffusion suppression coefficient — Kcff/'^Sp 
(where Kcs and Ksp are the effective and Spitzer diffu- 
sion coefficients) in the turbulent magnetized ICM. The 
obtained value turns out to be about 0.01 for a range of 
typical parameters of the ICM gas. This value is substan- 
tially lower than what have previously been assumed. We 
discuss implications of this result for the cluster properties, 
dynamics and formation. 

2. PARTICLE DIFFUSION IN THE CLUSTER HALO 

For simplicity we assume that the ICM gas is hydrogen 
and neglect the effect of helium. We also assume that the 
magnetic field in the ICM is highly turbulent and chaotic. 
Particle transport in such a system is strongly influenced 
by the field properties. 

A particle in the magnetic field gyrates about a certain 
field line and can move freely along it. Thus, on average, 
the particle closely follow the field line. A particle can 
jump through the distance ~ = vth_mpc/eB (the Lar- 
mor radius) to a neighboring field line in particle-particle 
collisions. Since the Larmor radius is many orders of mag- 
nitude smaller than the mean-free-path in the ICM, the 
cross-field diffusion is greatly suppressed. For future refer- 
ence, the Larmor radius, the proton mean-free-path for the 
typical conditions in the galaxy cluster halo, r ~ ^vir 1 

and 

the virial radius are tl = 2.3 x IQ-^^ j^jp^ j,i/2^ B^„q„ 

A = 1.3 Mpc T^^v A25\, rvir = 2.3 Mpc TgVev, where 
A = Up/up ~ 25 is the overdensity at the present 
epoch, z = 0, characterizing the ICM gas density at 
the virial radius, Up ~ 2 x 10~ ^ cm~^ is the mea n 
density in the Universe at z = (jSpergel. et al.l [20061 ) . 
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A25 = A/25, Tskcv = T/(5 keV), the ICM magnetic 
field i?inG = -8/(10^^ gauss). The thermal and Alfven 

speeds are vth = iikT/mpf/^ ^ 1.2 x 10^ T,V„i cm s"! 



•^SkcV 

VA = B/iATTmpUpY^^ = 9.7 x lO^BnoA^^^^ cm s'^. Be- 
low we consider three different scenarios for the particle 
transport: the case of static chaotic fields, the Alfvenic cas- 
cade and transport by large-scale turbulent motions (in- 
cluding the hydrodynam ic cascade]. We show that trans- 
port in Alfvenic cascade iLazarianl ( 2006f) is suppressed by 
some power of the Alfvenic Mach number, whereas the 
static fields and the fluid motions yield similar values for 
the diffusion coefficient for typical cluster conditions. 

Static chaotic fields — In a medium with strong mag- 
netic turbulence, the magnetic field lines themselves can 
be chaotic, such that the separation of two field lines 
increases exponentially with distance along them, d ^ 
exp(^||/ZB), with Ib being the field correlation length. 
The transport of a proton along such field lines re- 
sults in effective diffusion coefficient is suppressed com- 
pared to the classical Spitzer value Kgp — Awth — 4.8 x 

10'^^ cm^ s~^ ^5kcV ^25^- Note that here we shall not use 
the ambipolar rate, which is a factor of two largerQ 

If the magnetic field has a well defined, single scale Ib, 
the suppression factor due to th e field fine geometry, i.e., 
the field line tangling alone, is (jRechester fc RosenbluthI 
ll978HChandran fc Cowlev 19981) /field ~ [3ln{lB/rL)r^ - 
1/60. Here the factor 1/3 appears because the particle 
motion is effectively one-dimensional (i.e., along a field 
line), rather than three-dimensional. Here we also used 
the observat ional value of the typical fie ld correlation scale 
/s ~ 1 kpc (iVogt fc Enfilinl[2003l[200l . 

If Ib ranges through two decades or more in scale, as 
one can expect in the inertial range of strong magneto- 
hydrodynamic (M HD) turbulence, the suppress ion due to 
tangling is weaker ([Narayan fc Medvedevll200ll ): 



old 



1/5. 



(1) 
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In this case, the velisity on the injection scale is Vl 
For a general case Vl va, see iLazaria 

Another effect that is li miting the diffusion rate 
is the magnetic m irroring (Chandran fc Cowlev 'l998t 
iNaravan fc Medvedev 2001; Malyshkin fc Kulsrud 2001). 
Due to the conservation of the adiabatic invariant of a gy- 
rating particle, ^ = {mv'j_/2)/B, the particle may not be 
able to penetrate through the region of strong magnetic 
field, depending on the particle's pitch angle and the field 
strength. In clusters, the field inhomogeneities are large, 
SB ~ B, therefore the mirroring effect can be significant. 

Accurate calculation of the diffusion suppression due to 
mirroring in various regimes (single and multi-scale) and 
with various m odels of turbulence has been done in an ex- 
cellent work by (Malyshkin^KulsruS (IIMI). The mirror- 
ing suppression depends on the ratio X/Ib, an it is strong 
when Ib, as in the cluster halo. Note that the mirror- 
ing effect is small in t h e clus ter core where X < Ib, hence 
INaravan fc Medvedev! ([2001) argued that /mirroring <, 1- 
The results obtained by ^Mal vshkin fc KulsrudI ()2001h can 



be approximated quite well as follows: 

/mirroring =^niin{[ZB/(10'A)]'/', l}. (2) 

The overall diffusion suppression factor is 

/s = /field /mirroring ^ 4.3 X 10^"^. (3) 

Here we used ~ 1 kpc, and we assumed the most favor- 
able case of multi-scale turbulence with /field '^1/5. Note 
that even if the magnetic fields are correlated on a clus- 
ter size scale, 1 Mpc, the suppression is still rather 
strong: /s ~ 0.05. Thus, the magnetic fiel d suppressi on is 
much stronger than Jb ^ 1/3 assumed in (jLoebll2007t ). for 
which one have obtained that about 10% of the cluster gas 
will evaporate during the Hubble time. Simply re-scaling 
this result to the value of /b ~ 4 x 10"'^ yields the evap- 
orated fraction of order 0.1%. This result is essentially 
independent of the fraction of open field lines crossing rvir. 

Alfvenic turbulence — The particle diffusivity in the 
Alfvenic MHD cascade (b ut neglec t ing th e mirroring ef- 
fect) has been studied by LazarianI (|2006[ ). He has found 
that to sub-Alfvenic turbulence, Vgas/vA = Ma < 1, the 
diffusion coeffici ent is suppressed by a factor of M\ < 1 
compared to the INaravan fc MedvedevI (|200ll ) result: 



/f 



field 



(l/3)Ml, 



Ma < 1, 



(4) 



cf. Eq. (8) in 'LazarianI (|2006ll2007f ). The factor 1/3 above 
is approximate, howev er it is not very different fro m the 
factor 1/5 obtained bv INaravan fc MedvedevI (|200lD . 

For the super- Alfvenic turbulence, i.e., gas motions with 
the Alfven Mach number Ma > 1, there are two regimes. 
In the first one, the particle mean free path is large com- 
pared to the so-called "Alfven length", X > Ia, where 
I A ~ IbM7^ < Ib, and the suppression is 



/field il/3)ilA/X) ~ ilB/X)M^^, Ma > 1, (5) 

cf. Eq. (4) of iLazariaiil (|2006l ). and note that < A. In 
the second case, X < Ia < Ib, one recovers the /fioid '^1/3 
limit of one-dimensional diffusion. This case, however is 
not very relevant for cluster evaporation, because the mean 
free path is comparable to the system size and is much 
larger than the outer scale of turbulence, ~ 1 Mpc ^ Ib- 

As in any magnetic turbulence, the conservation of the 
magnetic moment results in additional suppression due to 
particle trapping, as discussed above, Eqs. (|2l3p . Since the 
"field" part of the suppression factor alone is diminished 
by some power of the Alfvenic Mach number in both sub- 
and super-Alfvenic regimes, we conclude that the static 
field model puts less stringent limits on thermal transport. 

Fluid turbulence — Turbulent fluid eddies can also con- 
tribute to the particle transport. The dynamica l diffusion 
coefficient has been also evaluated bv iLazarTa ^ (p006l : 



/^dyn 



Cd LVl, 
(3Cd LVlMI, 



Ma ^ 1; 
Ma < 1, 



(6) 



wher e /3 ~ 4 and Cd « 1/3 (|Lesieuijll990l ) (cf., it is not 2/3 
as in lLazariaiill2006l because we consider proton transport 
only). Note that transport in sub-Alfvenic turbulence is 



2 The ambipolar diffusion coefficient, Kamb = (Te -I- Tp)Ksp^ei^Sp,p/ (TpKSp,e + T'eKsp,p): where Ksp,e = Aefth.e and KSp,p = ^p^th.p < KSp,e 
for Te ~ Tp, following from plasma quasi-neutrality, is commonly used for two-component e~-p plasma. Here, we consider diffusion of a 
supra-thermal component of the proton distribution through gas, not the global escape of the plasma bulk, hence the Spitzer diffusivity is used. 
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suppressed by M\. Here L is the eddy scale and Vl is the 
turbulent velocity on this scale. 

It is instructive to represent the Spitzer diffusivity as 



KSp ~ 0.57 r5iiv^25^ (r-virWth), 



(7) 



i.e., it is very close to r-vir^th for typical cluster halo con- 
ditions. In the fluid turbulence, the combination LVl 
must be substantially smaller than rvir^th. Indeed, if 
LVl ^ Tvi-cVth, eddies of L fvir/2 and Vl ^ 2uth break 
the cluster apart: in virial equilibrium v^^^ = 2v1^ and gas 
parcels with Vl > are gravitationally unbound. 

Some fraction of gas can be lost in such violent events 
(e.g., mergers), but this cannot be called a steady-state 
evaporation. 

Numerical simulations of turbulence in clusters, tak- 
ing int o account decayin g turbulence, merger s and 
wakes (Dolag, et alj 120051 : iSubramanian. et al. I 120061 : 
IVazza. et al. 2006') yield a typical value of LVl ~ 4.5 x 
10" kpc km s-i ~ 0.016 T^^^ly (r^h-vth)- This yields 



Kdyn/KSp - 0.009 T^^^{^A25, 



(8) 



for Ma > 1, which is comparable to diffusion in static 
fields. Whether and how much the ICM fields will limit 
particle escape through Tvir has never been studied, but it 
goes beyond the diffusion approximation and may require 
a numerical approach. 

Diffusion length and evaporated gas fraction — Con- 
sidering several scenarios, we found that particle diffu- 
sion in the magnetized ICM is suppressed by at least a 
factor of a hundred or more, compared to the Spitzer 
value. Thus, the static field scenario represents an op- 
timistic estimate (lowest suppression) and we used it for 
further analysis. We now evaluate the diffusion length, 
through which a proton can travel during the Hubble time, 
iff = if-Mz=o - 4 X 1017 s: 



^diff = (^cffiif) 



1/2 



0.27 Mpcfl/'T^l^A-^/', 



(9) 



where /b,o.oo4 = /_b/0.004. Note that the diffusion co- 
efficient is inversely proportional to the local gas den- 



sity, which, in turn, is oc 



for the NFW pro- 



file (H avarro. et al.l Il997t ). Formally replacing A with 
A(r/rvir)~^, we see that the diffusion from the inner 
parts of the cluster is strongly suppressed. Hence, the 
highly coUisional gas at r < 0.5rvir cannot replenish 
the evaporated suprathermal component of the particle 
distribution. Because the escaping supra-thermal pro- 
tons move at speeds Vcsc ~ "n^'^vth ~ the 
diffusion distance above is accurate to a factor of or- 
der unity. Since most of plasma particles, both thermal 
and supra-thermal, are effectively trapped, their residence 
time in between magnetic mirrors can be comparable to 
the Hubble time. Supra-thermal particles with velocity 

V = ryi/^Wth can experience N collisions per the Hub- 
ble time: N ~ tu/TcoW ^ 12 (77 Tskov)"^^^ A25, where 
Tcoii = \{v)/v = 3.3 xlO^^ s (ryTgkev)^^^ A^/ is the particle 
collision time, ~ 2.3 for v be ing equal to the escape ve- 
locity, Vesc, at Tvir (|Loebll2007[ ) ■ Thus, there can be signif- 
icant re-distribution of trapped, v < Vcsc, and untrapped, 

V > Vcsc, components (i.e., an initially trapped supra- 
thermal proton may become untrapped and vice versa. 



due to collisions). This may change, perhaps lower, the 
effective 77. In this regime, numerical computation of the 
effective 77 using a self-consistent proton distribution ac- 
counting for diffusion, trapping and leaking of particles is 
highly desirable. 

The diffusion length is small fraction of the cluster size: 



^diff/fv 



n 19 f^^^ rp3/4 A —1/2 
"•^^ J_B, 0.004 -'skcV ^25 



(10) 



Thus, cluster evaporation effectively occurs from a thin 
layer of thickness ~ Idis O.lrvir. Particles from radii 
f < (''vir — ^diff ) cannot escape from the cluster within the 
Hubble time. 

The NFW density profile (jNavarro. et all I1997D is 

ngas(a;) = no/[x{l + x)^], where x = r/rg, = r^ir/c, c is 
the concentration parameter, uq is the normalization such 
that n{xvir) — fipA and Xvir = fvw/rs = c. The mass pro- 
file of the NFW cluster is M(a;) = Mo[ln{l + x)-x/{l + x)] 
with AIq = Airr'^nQmp. The fraction of the gas mass in the 
shell of thickness Zdiff to the total mass of gas in the entire 
cluster is 



hcU 



In 



c+l 

c+1-4 



(c+l)(c+l-5c) 



M, 



cluster 



ln(c+l)-c/(c-f 1) 



(11) 



The supra-thermal fraction of the proton distribution de- 
pends on the concentration parameter. It has been evalu- 



ated to be 0.064 < e < 0.092 at 



and for 2 < c < 8 



(|Loebll2007t ). Because of trapping of particles in magnetic 
traps, the particle distribution can evolve, therefore the 
actual value of e can differ. The fraction the evaporated 
gas is, therefore. 



AL 



hell 



6^CV(C+1)^ 



Mclustor ln(c-Kl)-c/(c+l) 



7.2 X 10" 



(12) 

where we used that ^ = 0.12,c = 4 and the corresponding 
e = 0.076. Thus, the evaporation fraction of the ICM gas 
is about 1%. Figures [1] and [2] represents the evaporated 
fraction /i as a function of the ICM temperature and the 
field correlation length, Ib, for several values of c and Ib- 

3. DISCUSSION 

We evaluated the suppression factor of the dif- 
fusive transport of protons in the ICM with mag- 
netic fields in three scenar ios: (i) static chaotic fields 
(|Naravan k Medvede"vll2001l) . (ii) Alfvenic and (in) fiuid 
turbulence ( LazariaiJ 20061 ). Note that magnetic fields are 
still required even in case (iii), otherwise the fiuid approx- 
imation breaks down, because A ~ Tvir. We considered 
diffusion suppression is due to (i) a purely geometric effect 
of turbulence and (ii) the dynamic effect of particle mirror- 
ing and trapping in inhomogeneous fields. We have found 
that particle diffusion is suppressed to at least 1% of the 
Spitzer value, and usually much more. For further analy- 
sis we used the suppression in static fields as an optimistic 
representative value. 

The smallness of the effective diffusivity, compared to 
the Spitzer value substantially lowers the overall evapo- 
ration gas fraction. We calculated the diffusion distance 
for typical cluster parameters and found that a particle 
can diffuse through the distance of about /diff ~ O.lrvir 
within the Hubble time. Thus, particles that are at radii 
r < [rvir — Idis) ~ 0.97'vir cau never pass through the virial 
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radius and leave the cluster. Consequently, cluster evap- 
oration occurs from the thin cluster "skin" of thickness 
~ O.lrvir ^ 200 kpc. The fraction of the gas evaporated 
from this skin can be rather high, perhaps about 10% o r 
even more, comparable to the estimate bv iLoebl ()2007t ). 
The exact value shall be obtained from self-consistent sim- 
ulations taking into account particle trapping, particle col- 
lisions, evolution of the distribution function and particle 
leaking out of the cluster. 

We make a prediction that the thermal energy deficit 
shall be the strongest at the outskirts of the cluster r ^ 
Tvir , in the shell of few hundred kpc in thickness. Deeper 
in the cluster, the thermal energy deficit shall be negligi- 
ble. Of course, at radii of few tens kpc, a cooling flow, if 
present, may lower the thermal energy as well. 

Since the diffusion distance is small compared to the 
cluster size, the fraction of the evaporated gas compared 
to the total cluster mass is small too. We evaluated it to 
be less than 1% for typical cluster conditions, T ~ 5 keV 
and ~ 1 kpc. Although the concentration parameter c 
can vary for different clusters, its effect on the evaporated 
fraction, is not very significant. The ICM gas tempera- 
ture and the field correlation scale do affect fi. Figures [101 
show how fi depends on T and / b ■ Interestingly, by boost- 
ing Ib to ^ 10 kpc and T to ^ 15 keV altogether, one can 
increase the evaporated fraction to about 2% only. One 
cannot exceed the value of ~ few % even by pushing I b 
to the maximum possible value of ~ 3 Mpc, — the cluster 



size. Apparently, the predicted evaporated gas fraction 
is below the observational value of 22 ± 10% by a large 
margin. 

Now, let's look at the result from a different point of 
view. Let's suppose that the observed thermal energy 
deficit of ~20% is due to cluster evaporation. What 
assumptions shall be relaxed in this case? We identify 
two alternatives. First, the ICM shal l be e ffectively un- 
magnetized, then the result of Loe3 ()2007f ) for the low- 
suppression scenario holds and the evaporated fraction can 
be of order 10%. In order to neglect the effects of mag- 
netic field, the particle Larmor radius shall be comparable 
to the size of a cluster, ^ rvh-. This sets the limit 
on the field strength: B < 10"^^ gauss. Second, if the 
ICM fields are present, they shall be relatively homoge- 
neous with SB <^ B. In this case, particle trapping is 
small and protons can easily escape along open field lines. 
The evaporated fraction will be reduced by the fraction of 
the cluster surface threaded by open fields lines. If this 
suppression is not strong, i.e., of order unity, then one can 
again expect ~10% value for the evaporated gas fraction. 
Alternatively, cluster mergers can lead to some gas loss. 
However, this process is episodic and is already taken into 
account in LSS hydro simulations. 

We thank the referee, Alex Lazarian, for discussion 
and comments. This work has been supported by grants 
NNG04GM41G, NNX07AJ50G and DEFG0204ER54790. 
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Fig. 2. — Fraction of the evaporated gas as a function of Ib, the maximum correlation length the magnetic field. Here it runs from the 
observed value of ~1 kpc to the maximum possible scale, the cluster size, ~3 Mpc. The solid curve is computed for T = 5 keV and c = 2, the 
dotted is for T = 15 keV and c = 2, and the dashed is for T = 5 keV and c = 8. 



